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ABSTRACT 

A  computer-controlled  system  for  automated,  remote 
measurement  of  phase  and  amplitude  shaded  doublet,  triplet, 
and  quadruplet  acoustic  arrays  beam  patterns  was  developed 
and  tested.  Experimental  data  were  collected  with  an 
anechoic  chamber  to  minimize  backscatter  and  other 
interference  effects. 

Preliminary  impedance,  reactance  and  resistance  plots 
were  made  on  several  speakers  to  ensure  that  a  matched  set 
with  similar  characteristics  was  utilized  for  the  arrays. 
Once  the  arrays  were  established,  a  computer  program  was 
developed  for  an  HP-86  computer  which  allowed  it  to  rotate 
the  arrays,  sample  and  collect  data,  and  plot  the  results. 

Various  combinations  of  amplitude  and  phase  shading  were 
applied  to  the  dipole,  tripole,  linear  quadrupole  and 
quadratic  quadrupole  array.  The  results  were  then  compared 
to  computer  generated  beam  patterns  which  agreed  to  i  2dB. 
This  uncertainty  was  mainly  generated  by  the  slight 
deviations  in  the  positioning  of  elements. 
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I.   INTRODUCTION 

A.  PURPOSE  OF  STUDY 

The  purpose  of  this  investigation  was  to  develop  and 
verify  a  computer-controlled  experimental  system  to  measure 
and  plot  the  pressure  fields  of  various  acoustic  arrays  in 
an  anechoic  environment.  This  was  to  be  done  using  various 
combinations  of  phase  shading  and  amplitude  shading  of 
doublets,  triplets,  and  quadruplets.  Experimentation  was 
conducted  under  computer  control  to  maximize  the  amount  of 
data  available  and  to  ensure  the  highest  degree  of  accuracy. 
Additionally,  a  comparison  was  conducted  with  the 
theoretical  predictions  obtained  by  Park  [Ref.  1], 

B.  BACKGROUND 

The  use  of  dipoles  and  quadrupoles  for  beam  forming  and 
direction  finding,  both  in  air  and  in  water,  has  been  a 
subject  of  intense  investigation  since  the  1940's.  The 
basic  radio  direction  finding  apparatus  of  the  1940's  was 
nothing  more  than  a  set  of  orthogonal  dipoles  which 
consisted  of  two  crossed  loops.  The  relative  signal 
amplitude  from  the  two  dipoles  was  used  to  determine  the 
direction  of  arrival  of  a  signal.   The  angle   is  equal  to 

Arctan  (y/x)  where  x  was  the  amplitude  of  one  dipole  and   y, 

c 
the   amplitude   of  the  other  dipole.    This  presented  a   180 
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bearing  ambiguity  which  was  resolved  by  comparing  the  phase 
of  the  X  and  y  dipole  signals  with  the  phase  of  a  co-located 
omni-directional  antenna.  [Ref.  2] 

In  the  1950's,  extensive  studies  were  conducted  using 
planar  arrays  with  multiple  acoustic  and  electromagnetic 
dipole  elements.  These  studies  led  to  the  formulation  of 
planar  arrays  employed  in  radar  systems  of  the  1960's. 
Further  research  resulted  in  the  ability  to  steer  beam 
patterns  by  use  of  multiple  elements,  phase  shading,  and 
amplitude  shading. 

In  1976,  work  was  done  by  Hughes  and  Thompson  to 
investigate  acoustic  communications,  particularly  in 
underwater  acoustics  for  cases  where  the  position  of  the 
source  or  the  receiver  is  fixed.  They  investigated  the 
tilted  directional  response  patterns  below  the  horizontal 
plane  due  to  amplitude  weighting  and  a  single  90  phase 
shift.  Here,  a  method  of  achieving  a  tilted  pattern  from  a 
symmetrically  shaped  array  without  using  numerous  phase 
shifts  or  time  delays  was  discussed.  [Ref.  3]  Moses  and 
Smith  also  examined  beamforming  this  time  with  a  cylindrical 
array  using  real  (rather  than  complex)  shading  coefficients 
[Ref.  4]. 

In   1985,   Flanagan  examined  the  use  of  computer-steered 

microphone   arrays  which  were  two-dimensional  planar   arrays 

of  9  X  7  elements.    Here  an  array  was  designed  and  computer 

steered  to  discriminate  against  sound  arrivals  from  all   but 

the  desired  directions.  [Ref.  5] 
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One  example  of  the  Navy's  current  use  of  dipole 
technology  is  in  sonobuoys .  The  SSQ-53  DIFAR  Sonobuoy 
possesses  both  an  omnidirectional  and  bi-directional 
hydrophone  capability.  The  SSQ-77A  VLAD  Sonobuoy  possesses 
an  upper  and  lower  array  strung  in  a  single  line  array. 
Beamforming  is  accomplished  by  the  use  of  both  phase  shading 
and  amplitude  shading. 

These  advancements  in  the  study  of  arrays  ranging  from  a 
single  dipole  to  large  scale  planar  arrays,  coupled  with 
computer  control,  form  the  basis  for  this  further  study  of 
shaded  arrays.  There  is  currently  a  related  research  effort 
ongoing  in  the  study  of  computer  controlled  arrays  and  beam 
patterns  being  conducted  by  the  Environmental  Physics  group 
at  the  Naval  Postgraduate  School. 

C.   REPORT  OVERVIEW 

In  Chapter  II,  a  basic  introduction  to  theoretical 
pressure  and  directivity  calculations  of  a  dipole  is 
developed.  This  is  then  used  as  a  base  to  expand  to  the 
physical  case  and  the  electrical  linear  systems  case. 
Finally,  these  results  are  applied  to  different  forms  of 
quadruplets  and  triplets  to  show  theoretically,  the 
expected  results. 

Chapter  III  contains  a  summary  of  the  beam  pattern 
measuring  system,  the  computer  program  development, 
experimental  procedure,   and  representative  data  and  plots. 


13 


Additionally,  it  provides  a  description  of  the  system  which 
was  assembled  for  measuring  the  pressure  fields  of  the 
various  arrays. 

Chapter  IV  contains  a  summary  of  data  analysis  and 
verification.  Additionally,  representative  plots  are 
provided  for  amplification  of  arguments. 

Chapter  V  contains  concluding  remarks  and 
recommendations  for  further  work. 

Appendix  A  provides  a  brief  description  of  each  piece  of 
equipment  and  summary  of  its  capabilities  as  applied  to  this 
experimental  investigation. 

Appendix  B  lists  the  BASIC  computer  language  program 
used  for  data  acquisition  and  formulation  in  this  study. 

Appendix  C  contains  specific  plots,  both  experimental 
and  theoretical,  obtained  as  a  result  of  this  investigation 
and  that  of  Park  [Ref.  1].  These  plots  were  obtained  for 
beam  patterns  of  various  array  configurations. 
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II.   THEORY 

The  dipole  has  been  extensively  explored  and  is  now 
utilized  as  a  basic  unit  in  the  construction  of  a  great  many 
arrays  and  antennae.  In  this  section  we  will  elucidate  the 
theory  using  two  approaches  1)  an  acoustic  approach  in  2 
dimensions  and  2)  a  linear  systems  approach  in  three 
dimensions.  These  two  approaches  are  related  to  each  other 
and  yield  the  same  result  after  minor  changes  are  made  in 
the  geometry  of  the  problem.  Both  are  presented  because  of 
direct  Navy  applications  dealing  deal  with  SONAR  arrays. 

Following  the  presentation  of  each  of  these,  the 
equations  for  the  triplet  and  quadruplet  arrays  utilized  in 
this  research  effort  will  be  shown. 

A.   DIPOLE  RADIATION-ACOUSTIC  APPROACH 

Because  it  is  usually  desirable  and  useful  to  be  able  to 
determine  the  direction  from  which  an  acoustic  signal 
originated,  arrays  were  developed  to  provide  this 
information;  a  simple  example  of  this  is  the  acoustic 
dipole.  A  simple  dipole  consists  of  two  simple  sources  of 
equal  strength  but  opposite  phase,  separated  by  a  distance  d 
Fig.  (2.1). 
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P(r9. ») 


Figure  2.1   Geometry  used  in  Deriving  the 

Radiation  Characteristics  of  an 
Acoustic  Dipole. 


The  pressure,   p(r,Grt),   measured  at  the  field  point  (r,0) 

due  to  source  1  and  2  respectively  are 

Pj  =  [A;(r  +  Ar,)]  exp[j((Ot  -  k(r  +  Ar,)]  ^j 

P^  =  (-A'(r  -  Ar,)]  exrfj{(ot  -  k(r  +  Ar,)]  (2) 

where  A  is  the  pressure  amplitude  of  an  individual  source 
measured  at  one  meter  from  the  source,  r  is  the  distance 
from  the  field  point  to  the  midpoint  between  the  sources, 
0)— 2  77 f  is  the  angular  frequency  and  k  is  the  wave   number. 

The  minus  sign  associated  with  P2  accounts  for  the  180 
phase  difference.   The  total  acoustic  pressure  is  the  sum  of 
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Pl  and  PX' 

P(r.e.t)r(A  r){{[exr(-jkArj)][lt(Ar^,'r)])-{[exp(-jkAr2)]'[l»(Ar,/r)]))exp(j(tot-  kr))   (^^ 
Considering  only  field  measurements  where  the  separation  d  is 


small  compared  to  the  distance  r,  we  can  use  the  far  field 
approximations,   ^^1,2  ^  (d  2)  sin  0      and     Ar.  ^/r  <  <  1 
This  yields  the  acoustic  dipole  radiation  pressure  field, 

P(r,e,t)  =  -j  (2A;r)  sin  ((kd;2)sin  0)}  exp  i(0)t  -  kr)  ^4^ 


180  e 


Figure  2.2  Cartesian  and  Polar  Representation 
of  Pressure  Amplitude  Distribution 
of  Dipole  Radiation.  [Ref.  6] 
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As  noted  in  Eq.  (4)  and  seen  in  Fig.  (2.2),  plotting  the  far 
field  for  d<<Af   pressure  nulls  exist  when  (kd/2)  sinQ  =77n 
where  n  is  an  integer.  [Ref.  6] 

The  twin  lobes  and  nulls  of  a  dipole  result  in  an  180 
ambiguity  when  trying  to  localize  a  source.  If  the  dipole 
is  combined  with  an  omnidirectional  receiver,  the  ambiguity 
can  be  resolved.  In  Fig.  (2.3),  various  combinations  of 
dipole  and  omnidirectional  receivers  are  combined  to  provide 
a  cardioid  response. 


BlOinCCTlOMAL 


NONCMOCCTIOMIU. 


COMBINATION 


9-cD4^ 


IZO        0*-<l 


Figure  2.3  Directional  Diagrams  of  Various  Combinations 
of  Bidirectional  and  Nondirect ional 
Microphones.  [Ref.  7] 


Here  R,    is  the  voltage  output  of   the  omnidirectional 
receiver  and  R2  the  voltage  output  of  the  directional  dipole 
receiver.    It  is  not  important  for  the  two  sources  to  be  of 
the  same  sensitivity,  but  the  ratio  R^/Rx  must  be  unity  for 
true  cardioid  response.    The  phase  shift  between  these   two 
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receivers,  which  yields  directionality,  is  achieved  through 
a  difference  in  the  phase  characteristics  of  the  omni  which 
is  differentiated  (high  pass  filtered)  to  compensate  for  the 
linear  frequency  dependence  of  the  dipole  sensitivity.  [Ref. 
7]. 

This  design  is  presently  used  in  SSQ-77A  VLAD  Sonobuoys 
where  a  cardioid  is  formed  pointing  in  the  downward 
direction.  This  maintains  the  null  in  the  upward  direction 
in  order  to  alleviate  high  frequency  interference  from  wind 
and  sea  surface  noise. 

B.   LINEAR  SYSTEMS 

Consider   the  geometry   for  a  dipole  as  shown   in   Fig. 
(2.4)  : 


ip.9,i) 


Figure  2.4  Dipole  Geometry  [Ref.  8] 
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In  this  representation  as  developed  in  Ref.  (7),  C^  and  C_i 
are  complex  weighted  point  sources  which  are  free  to  rotate 
in  the  x-y  plane  and  d  is  the  distance  between  elements. 

The  complex  weighted  response  factor  ci  ^^  ^^^   ^^^  point 

source  is 

S  =  ^i  exp(j(p.)  ^^^ 

where  ai  is  the  amplitude  and  ^i  is  the  phase  angle.  For  a 
dipole,   the  complex  weighted  response  factors  are   Cj  =  a. 

and  '^-1  =  -a^  • 

For   a  wave  originating   from  a  distant   source  and 

propagating  to  the  receiver,  it  is  observed  that  there  is  an 

essentially  constant   pressure  magnitude  over   the   dipole 

aperture.   Because  the  phase  varies  with  position  along  the 

wave  front  in  a  traveling  wave,   the  response  of  the  dipole 

is  given  as  a  function  of  frequency  and  direction  cosine  u 

where  u  =  sin0cos^.  The  response  is 

D(r,u)  =  c,  exp{j7rdu;>.)  +  C2  exp(-i7idu;X)  (6) 

This  equation  is  simplified  for  the  dipole  to: 

D(f,u)  =  2iaj  sin{7rdu;X)  (7) 

Since  u   is  a  direction  cosine,   it  is  only  valid  over  the 

interval  -l<u<l.   As  a  result,  0<G<7r   and  ()^\\f^2K   .  u  =  o 

implies   a  broadside  or  on-axis  geometry  and  u  =   1   implies 

an  end-fire  geometry.   A  plot  of  D(f,u)   as  a  function  of  u 

can  be  used  to  approximate  the  response  of  a  dipole.  [Ref. 
8] 
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An  example  of  this  two  element  response  of  two  point 
sources  with  a  separation  d  equal  to  one  half  the  wavelength 
of  the  traveling  wave  results  in  a  real  response  for  values 
of  u  from  (->^;2d)  to  ( +  X;2d)  •  If  the  sources  are  assigned 
dipole  weightings,  substitution  of  the  f unctions  G  and  N^*  for 
direction  cosine  u  result  in  the  unnormal ized ,  far  field 
beam  pattern 

D(r,0.vi/)  =  i2a|  sinl(7rd;X)  sin  6  cos  \\f]  (8) 

This  linear  systems  theory  can  be  easily  applied  to  the 
axial  quadrupole  and  Tesseral  quadrupole.  The  axial 
quadrupole   is  an  array  of  three  complex  weighted  point 

sources  where  c,  =  a,  =  c.j  =  a.,  =  1     and  0^  =  3^=  -2  as  in 
Fig.  (2.5). 

Here  C  can  also  be  thought  of  as  2  sources  with  each  having 
an  amplitude  weight  equal  to  -1. 


ip.8.f) 


Figure  2.5   Axial  Quadrupole  Geometry  [Ref.  8] 
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with        the      1,-2,1      amplitude      weighting,        the        far-field, 
unnormalized   directivity   function    is   defined   as: 

D(r.r^)  =  [exp(-j2;Tf^d'2)l  -  2  4-  [exp(j27rf;d2)l  (9) 

where  f^  =  u  X  =  sinGcosvi/  X  .      This    reduces   to: 

D(r.G.H/)  =  2  [cos{(7rd;>.)  sinGcosH/)  -  1]  (10) 

In  the  case  of  a  Tesseral  quadrupole  Fig.  (2.6),  a 
planar  array  of  four  amplitude-weighted  point  sources,  the 
unnormalized  directivity  function  is: 

D(r,r^.r  )  =  2Q  cos[Trd(f^  +  r,)]  -  2Q  cos[7rd(f;  -  f,)] 

where      f^  =  u  X  =  sinGcosij/';^  r      f   =  v/T.  =  sinGsinV|/  X 

direction  cosine  v   =   sin  0  sin  t .   This   yields 

D(r,G,y)  =  -4Q  sin[(7rd  X)  sinGcosy]  sin[(7rd;>.)  sinGsinxj/]  (12) 


and 


(11) 
the 


ip.S.f) 


Figure  2.6  Tesseral  Quadrupole  Geometry 
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which  is  the  far  field,  unnormalized  directivity  function 
for  this  quadrupole  as  a  function  of©  and  r .  [Ref.  8]  The 
above  development  is  particularly  important  because  this 
type  array  was  investigated  in  the  course  of  this  work. 

C.   GENERAL  QUADRUPOLE 

Park,  in  his  thesis  (Ref.  1),  developed  an  equation  for 
calculation  of  a  four  element  array  Fig.  (2.7).  The 
equation  was  developed  by  adding  complex  pressures  from  four 

sources  to  yield 

P  =  _L  eJ(wt  -  kR  -  kDjSinGsind))  ^  _^  Jim  -  kR  +  kD^sinesinO)  -  (p,)      /,  ox 
R  R  -       V      ; 

+  JL  ei("^  -  ^^  +  kD3sinecos(I>  -  ip^) 
R 

4.     ^^    gj(wt  -  kR  -  kD^^sinGcosO)  -  (p^) 
R 

where  [( 1  4- A4- B-+-C)  R]  explj(cot  -  kR)]  is  the  maximum  possible 
pressure  amplitude  and  the  distances  D^^  D2  D^  D^  and 
measured  from  the  origins  to  the  particular  element.  The 
directivity  is: 

11(0,0))  =  {[cos(kD,sinGsin<I>)  +  Acos(kD,sinGsinCl)  -  (pJ  ^^^^ 

(1+A  +  B  +  C)  ^  '  ■^ 

+  BcosikDjSinGcosO)  -  (P3)  +  Ccos(kD_jSinOcos(I>  +  9^)]- 

+  lsin(kD,sinGsiii<I))  -  Asin(kD,sinGsin<D  -  (p,) 

-  Bsin(kD3siiiOcostI>  -ip^)  +  CsinlkD^sinGtosO)  +  (p4)l-)  "" 
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Figure  2.7   Quadrupole  Geometry 

If  D-^  =  D2  =  D3  =  D4  then  comparison  of  Figs.  (2.7)  and 
(2.6)  shows  that  the  configuration  is  simply  the  Tesseral 
quadrupole  of  Fig.  (2.6)  where  ^has  been  shifted  45  . 

Additionally,  the  axial  quadrupole  of  Fig.  (2.5)  is 
obtained  using  Park's  program  by  placing  source  B  at  an 
inf initesimally  small  distance  from  the  array  center 
(i.e.,  D3  =  0.00001  meters  for  any  A  )  with  amplitude 
weighting  equal  to  twice  sources  1  and  A  and  source  C 
amplitude   weighted  to  zero,   which  removes  its  contribution 
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from  the  array.  This  axial  quadrupole  can  also  be  obtained 
by  placing  sources  B  and  C  at  the  array  center  and  amplitude 
weighting  equally. 

These  equations  were  utilized  to  provide  verification 
data  for  the  laboratory  experiments.  The  three  element 
array  data  were  verified  by  amplitude  weighting  element  C  of 
Eq.  (14)  to  zero.  Dipole  data  were  verified  by  amplitude 
weighting  elements  B  and  C  of   Eq.  (14)  to  zero. 
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III.   EXPERIMENTAL  PROCEDURE  AND  DEVELOPMENT 

A.   BEAM  PATTERN  MEASURING  SYSTEM  OVERVIEW 

The  beam  pattern  measuring  system  is  an  instrumentation 
package  which  was  designed  to  power  and  rotate  an  acoustic 
array  and  then  receive  the  signal,  compute  the  pressure 
distribution  of  the  array  and  display  this  data  in  both  a 
tabulated  format  and  plot  of  sound  pressure  level  in 
decibels  versus  degrees  of  rotation.  All  data  were  taken  in 
an  anechoic  chamber  to  minimize  surface  interference. 

The  system  used  for  measuring  beam  patterns  of  different 
arrays  is  divided  into  two  major  groups.  These  are  the 
signal  generation  and  input  section  and  the  automated 
advancement,  data  collection  and  display  section.  An 
overview  of  these  is  now  presented  along  with  a  block 
diagram  of  the  system.  Fig,  (3.1).  A  detailed  description 
of  each  piece  of  equipment  is  contained  in  Appendix  A. 

The  Hewlett-Packard  (HP)  Function  Generator  Model  3314A 
was  utilized  for  all  signal  generation.  The  signal  passed 
through  an  individual  HP-467A  Power  Amplifier  to  each  of  the 
University  Sound  Model  ID-30  acoustic  driver  units.  Setting 
the  output  levels  of  the  driver  units  was  accomplished 
utilizing  the  variable  gain  control  on  the  power  amplifiers. 
This  constituted  the  signal  generation  and  input  section. 
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Figure  3.1  System  Diagram 
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The  receiving  and  processing  section  was  controlled  by 
the  HP-86B  computer.  It  directly  controlled  the  array 
rotation  and  the  reading  of  all  measurements  of  peripheral 
instruments  via  the  Hewlett-Packard  Interface  Bus  (HPIB) . 
The  computer  sent  clockwise  stepping  commands  to  the  motor- 
power/direction  controller  via  the  HP-3421A  Data 
Acquisition/Control  (DAC)  to  rotate  the  array  in  stepped 
increments  of  about  one  degree.  Following  each  advancement 
of  the  array,  the  computer  read  both  the  resistance  output 
of  the  array  rotation  measurement  potentiometer  and  the 
voltage  output  of  the  General  Radio  Model  1962  Microphone. 
The  array  rotation  potentiometer  was  a  0  to  10  Kohm 
potentiometer  mounted  on  top  of  the  rotation  motor  and 
measured  the  angle  of  rotation  of  the  array  (0  to  10  Kohms 
corresponded  to  0  to  360  ) .  These  resistance  and  voltage 
readings  were  then  stored  in  the  computer  for  use  in  the 
data  analysis  and  plotting  routines.  Once  the  data  was 
taken  for  360  of  rotation,  it  was  used  to  compute  the 
rotation  angle  and  its  corresponding  decibel  (dB)  level  for 
each  datum.  These  values  were  then  output  to  the  computer 
monitor  and  plotted  in  the  form  of  a  polar  plot  of  dB  versus 
degrees  of  rotation. 

B.   BEAM  PATTERN  PROGRAM 

The  Beam  Pattern  program  listed  in  Appendix  B,  consists 
of   a   series   of   sub   routines.    The  program  starts   by 
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providing  addresses  for  the  two  multimeters  and  the  data 
acquisition/control  unit.  Next  the  initialization 
subroutine  resets  and  initializes  these  three  instruments. 

Once  the  equipment  is  ready  to  operate,  the  incrementing 
subroutine  is  called.  This  is  an  iterative  and  interactive 
routine  used  to  rotate  the  array  and  measure  the  voltage  of 
the  sensed  pressure  field  and  the  corresponding  resistance 
reading  which  is  proportional  to  the  angular  position  of  the 
array.  Up  to  360  data  points  can  be  collected  per  data  run. 
After  the  initial  datum  is  read,  the  array  rotator  motor  is 
energized  for  0.5  milliseconds  to  increment  the  array 
azimuth.  The  array  is  then  allowed  to  coast  down  for  4 
seconds  prior  to  reading  the  voltage  and  resistance.  This 
selection  of  rotate  and  coast  down  times  is  discussed  later 
in  the  chapter.  Once  the  resistance  and  voltage  readings 
are  taken,  they  are  stored  in  an  array  for  later  transfer  to 
a  data  disk  and  for  use  in  the  plot  subroutine. 

Finally  the  plot  subroutine  calculates  the  relative 
decibel  level  of  the  beam  pattern  based  on  the  measured 
voltages.  It  then  takes  the  resistance  reading  and  converts 
it  to  the  angle  of  rotation.  The  calculated  decibel  levels 
are  then  normalized  to  zero  decibels  based  on  the  highest 
received  voltage  reading.  Finally,  the  data  points  are 
plotted  on  a  polar  plot  of  relative  sound  pressure  level  in 
decibel  versus  angle  of  rotation. 
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All  raw  data  is  stored  on  a  data  disk  after  a  complete 
data  run.  Each  data  run  is  stored  by  run  number,  as  shown 
in  Table  II.  It  can  be  retrieved  and  printed  at  the 
printer  as  received  voltage  and  resistance  for  each  datum 
point  or  printed  and  plotted  as  decibel  versus  angle  of 
rotation. 

The  physical  weight  of  an  array  was  a  major  variable  in 
this  program  with  regard  to  the  selection  of  the  number  of 
data  points  needed.  Arrays  weighing  less  than  4  pounds 
required  360  data  points  for  a  full  360  of  rotation.  The 
number  of  data  points  required  for  larger  arrays  was 
decreased  as  array  weight  increased.  For  the  four  element 
array,  weight  was  substantial  enough  (19  pounds)  to  require 
only  317  data  points  for  a  full  360  rotation.  This  was 
caused  by  the  fact  that  the  heavier  arrays  would  coast-down 
for  a  larger  distance  after  the  0.5  millisecond  motor-on 
period. 

Another  contributor  to  variation  in  the  required  number 
of  data  points  was  array  balance.  When  utilizing  a  three 
element  array,  the  array  was  sufficiently  out  of  balance  to 
require  a  greater  number  of  data  points  for  full  rotation. 
This  was  attributed  to  minor  binding  of  the  rotator  shaft  as 
the  array  coasted  down. 

Selection  of  motor  run  time  and  coast-down  time  was 
based  on  experimental  investigation.  A  0.5  millisecond 
motor  on  time  provided  1   rotation  regardless  of  coast-down 


30 


time  with  a  light  weight  array.  The  4  second  coast-down  was 
selected  because  smaller  periods  did  not  allow  the  array  to 
stop  and  the  resistance  multimeter  to  stabilize  prior  to  it 
being  read.  Motion  at  the  variable  potentiometer  when  the 
resistance  values  were  read  yielded  erratic  reading  on  the 
multimeter  due  to  electromagnetic  interference  generated  by 
the  motor, 

C.   EXPERIMENTAL  APPARATUS  DEVELOPMENT 

The  experimental  apparatus  was  designed  to  optimize  a 
thorough  search  of  two,  three  and  four  element  arrays  and  to 
establish  their  respective  beam  patterns.  Previously, 
several  experimental  investigations  have  been  conducted 
utilizing  a  two  microphone  array  with  a  single  source.  It 
was,  therefore,  decided  that  reciprocity  would  be  invoked 
and  the  array  would  consist  of  multiple  ID-30  driver  units 
(Appendix  A)  with  a  single  microphone  utilized  for  sensing 
the  pressure  field. 

Six  ID-30  driver  units  were  checked  individually  in  the 
anechoic  chamber  for  their  beam  pattern  impedence, 
resistance  and  reactance  characteristics.  This  data  was 
utilized  to  select  the  best  matched  set  of  four  driver 
units.  All  drivers  were  examined  between  5  Hz  and  3000  Hz 
in  an  effort  to  locate  and  identify  the  resonances  in  the 
individual  speakers.  Fig.  (3,2)  shows  a  representative 
resistance   curve  for   these  drivers  between  the  above 
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frequencies.  It  was  desired  to  conduct  all  analysis  at 
frequencies  in  relatively  flat  areas  away  from  resonances. 
This  was  done  to  assure  phase  stability  since  reactance  is 
not  varying.  The  speakers  utilized  in  this  investigation 
had  resonances  as  per  Table  I. 

TABLE  I 
ID-30  DRIVER  UNIT  RESONANT  FREQUENCIES 
Speaker  Number  ;  First  Resonance  (Hz) ;  Second  Resonance  (Hz) 

1  689  2290 

2  660  2225 

3  780  2150 

4  798  2129 

5  740  2253 

6  778  2125 

Matched  set  number  1  consisted  of  speakers  1  and  2  and  set 
number  2  of  speakers  3  and  4.  All  array  beam  patterns  were 
run  at  300  Hz  and  960  Hz  which  fell  in  the  flat  areas  of  the 
impedance  curves. 

The  array  existed  in  two  individual  configurations,  a 
linear  and  a  crossed  axis  system  with  variable  tilt  (Fig. 
(3.3)).  The  linear  system  was  used  for  investigation  of 
doublets,  some  triplets  and  quadruplets  (in  the  form  of  an 
axial  quadrupole) .  The  crossed  axis  array  proved  to  be 
extremely  versatile  and  useful,  providing  four  arms,  each 
70   centimeters   in   length,    for   holding   the   drivers. 
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Additionally,  it  provided  a  means  to  vary  the  angle©  a  full 
90°.  Each  driver  was  equipped  with  a  harness  assembly  which 
mounted  to  the  array  rod  in  such  a  manner  that  the  driver 
face  was  located  2.5  centimeters  off  the  axis.  This  did  not 
prove  to  be  a  matter  of  concern  since  all  pressure 
measurements  were  made  in  the  far  field. 


^/»©,y; 


Quadrtiplet 


Figure  3.3  Array  Configurations 
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A  scattering  and  shading  problem  of  small  magnitude  was 
encountered  because  of  the  large  center  piece  on  the  crossed 
axis  array.  This  scattering  problem  was  largely  resolved  by 
wrapping  the  center  piece  with  fiberglass. 

The  General  Radio  Model  1962  microphone  was  located  3.2 
meters  from  the  center  of  the  array  placing  it  in  the  far 
field.  It  was  mounted  at  the  same  height  in  the  anechoic 
chamber  as  the  array  center.  Changes  in  tilt  angle  theta 
were  accomplished  by  tilting  the  array. 

The  outputs  of  the  ID-30  driver  units  were  balanced 
prior  to  each  data  run.  It  was  done  by  aligning  the  array 
in  the  vertical  plane  with  the  drivers  facing  the 
microphone.  Then  the  output  of  each  driver  was  adjusted 
utilizing  the  variable  gain  control  on  the  power  amplifier. 
This  allowed  all  the  drivers  to  be  driven  within  a  received 
voltage  of  0.002  VAC  of  each  other  out  of  a  typical  level  of 
1  Volt. 

Output  signals  from  the  HP-467A  Power  Amplifiers  were 
monitored  on  an  oscilloscope  as  was  the  received  signal  from 
the  microphone.  Because  the  signals  were  clean  and  without 
any  substantial  noise,  filters  were  not  utilized  in  the 
microphone  circuit. 
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IV.   DATA  ANALYSIS  AND  VERIFICATION 

After  establishing  the  ability  to  accurately  rotate  an 
acoustic  array  while  collecting  and  displaying  sensed 
pressure  levels  under  computer  control,  it  was  necessary  to 
ensure  the  data  were  accurate.  This  could  be  done  by  taking 
a  set  of  parameters  and  obtaining  both  theoretical  and 
experimental  beam  patterns  for  them.  Verification  was 
conducted  utilizing  the  results  of  a  FORTRAN  program 
developed  by  Park  (Ref.  1)  which  yields  a  variety  of  useful 
plots  with  the  help  of  DISSPLA. 

A.   THEORETICAL  BEAM  PATTERN  PREDICTION 

With  the  program  developed  by  Park,  quick  beam  pattern 
plots  could  be  made  to  compare  with  experimental  patterns. 
As  an  example.  Fig.  (4.1)  depicts  a  three  dimensional 
representation  of  a  dipole.  Clearly,  this  particular  is  not 
of  much  use. 

The  three  dimensional  plot  was  sectored  for  greater 
utility.  Fig.  (4.2)  shows  a  sector  representation  for  a 
four  element  array  with  kDj  =  kD2  =  3.52  and  ko^  =  kD^  =  0. 
Note  that  the  kO-j^  utilized  in  the  remainder  of  the  thesis 
related  directly  to  the  spacing  from  the  center  of  the  array 
to  the  element  as  related  to  Fig.  (2.7).  Refer  to  Fig. 
(2.7).    The  Z-axis  of  the  plot  is  beam  pattern  (dB) ,  the  Y- 
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DIRECTIVITY  PATTERN 


kDl=kD2=  3.52 
kD3=kD4=  0 

D1=0J2  D2=02 
D3=0.   D4=0. 
VIEW  ANGl£=  20.  20 
P2=180  P3=0  P4=  0 
A=l.  B=0.  C=0. 


Figure  4.1     Three  Dimensional  Dipole  Representation 
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kD2=kD2=kD3=kD^=  10.55 

Dl=0.6  D2=0  6 
03=0.6  D4=0.^ 
VIEW  ANGLE=-h20.+60 
P2=  0  P3=  0  P4=  0 
A=l.  B=l.  C=l. 


Figure  4.2  Sectored  Three  Dimensional  Plot  (  4  Element  ) 
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axis  psi  (degrees) ,  and  the  X-axis  theta  (degrees) . 
Directivity  is  scaled  from  20  to  -40  dB  (with  0  dB  being  the 
highest  achievable  value  because  of  normalization) ,  psi  from 
180  to  360  ,  and  theta  from  0  to  -90  .  This  plot  appears 
to  be  very  busy  but  proves  to  be  quite  useful.  It  is  easy 
to  evaluate  this  sectioned  representation  by  following  the 
contour  line  of  the  directivity  above  the  line  theta  equals 
0  for  psi  from  180  to  360  .  This  provides  the  same  dB 
levels  at  180  <  psi  <  360  as  the  polar  plot  representation 
plotted  in   Fig.  (4.3). 

B.   EXPERIMENTAL  DATA-SINGLE  ELEMENT  DIRECTIONALITY 

As  stated  in  Chapter  III,  testing  was  conducted  at  300 
Hz  and  960  Hz  which  were  in  the  flat  areas  of  the  drivers' 
resistance  curves,  away  from  their  resonance  peaks.  This 
was  done  to  assure  phase  stability  between  drivers  since  the 
reactance  is  not  varying  (large  phase  shifts  do  not  occur 
except  near  resonance) . 

Initial  data  runs  were  conducted  with  a  single  ID-30 
acoustic  driver  to  establish  the  directionality  of  each  of 
the  drivers.  The  drivers  were  rotated  360  at  three 
separate  tilt  angles  of  the  driver  face,  0=0  ,  45  and  90. 
A  representative  graph  of  the  directionality  is  contained 
in  Fig.  (4.4)  which  was  made  with  the  driver  at  a  tilt  angle 
of  0  and  pointing  at  the  microphone  prior  to  rotation.  The 
dotted  line  is  the  laboratory  beam  pattern.   Note  that  there 


39 


PHI  VS.  DIRECTIVITY 
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THCTA=  0  DEORTT^ 

kDi=kD2=k03-kC4=  10.55 

1)1=0.6  D2=0.6 
D3=0.6  D4=0  6 
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Figure  4.3  Polar  Plot  of  a  Four  Element  Array 


40 


m 

'TJ 


Li.] 
h- 


III 
en 

LiJ 

! 

o' 

n. 


V 


/^ 


...  -O- 


X 


/ 


X 


/ 


/ 


\ 


\ 


h     43-    ft    -^ ]- 

I'  II  11/ 


\ 


— i- 


/  / 


X 


\.. 


\ 


-I- 


^ 


■/ 


■V 


•y 


>, 

■p 

•H 

1— 1 

(U 

C 

o 

•H 

\ 

+J 

■\ 

N 

u 

\ 

X 

o 

■H 

'.  ^ 

CO 

Q 

en 

n 

C3 

a 
a 

or. 
u. 

m 

> 

■H 

>^ 

Q 

:  1 

n 

1 

U 

.1 

7 

a 

-H 

to 
a 

. 

0 

,*» 

on 

u 

•'• 

II 

< 

/ 

rt 

o 
n 

1 
Q 

M 
U 

a 
en 

41 


is  a  1  to  2  dB  drop  below  the  0  dB  level  as  the  driver 
points  away  from  the  receiver  which  shows  the  driver  is  not 
completely  omnidirectional.  It  should  also  be  noted  that 
the  program's  normalization  routine  sets  the  highest 
received  microphone  voltage  level  to  0  dB  and  normalizes  all 
other  voltage  levels  in  the  beam  pattern  to  this  level. 
This  was  done  to  provide  greater  ease  in  plotting  of  beam 
patterns  so  as  to  produce  a  uniformity  between  plots. 

C.   EXPERIMENTAL  AND  THEORETICAL  DATA  -  VARIOUS  ARRAY 
CONFIGURATIONS 

An  investigation  of  two,  three  and  four  element  arrays 
was  conducted.  For  each  of  these  configurations,  several 
different  combinations  of  spacing  and  amplitude  and  phase 
shading  were  looked  at.  Table  2  at  the  end  of  this  chapter, 
is  a  list  of  these  variations  for  which  beam  patterns  were 
obtained.  Note  that  phase  shading  was  limited  to  either  0 
or  180  .  In  viewing  the  figures,  the  term  "x-pole"  is  used 
loosely  to  represent  the  term  "x-let"  (i.e.,  quadrupole  for 
quadruplet) . 

For  the  two  element  array,  beam  patterns  were  run  in  the 
laboratory  for  both- in- phase  and  180°  out-of-phase  doublets. 
An  example  of  the  out-of-phase  doublet  beam  pattern  obtained 
in  the  anechoic  chamber  is  shown  in  Fig.  (4.5).  (Note  that 
in  this  figure  kd  is  the  quantity  2kD|  when  Dj  =02-) 
Its  theoretical  verification,  obtained  from  Park's  program 
is   shown  in  Fig.   (4.6).   A  contributing  factor   to  their 
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THETA  VS.  DIRECTIVITY 

0 


PHI=  90  DEGREE 

kDi=kD2=  5.276    ,      kD3=kD4=  0 

D1=D2=  0.5 

D3=D4=  0 
P2=180  P3-  0   P4=  0 
A-1.   D-0.   C=0. 


Figure  4.6  Theoretical  Dipole  Plot 
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small  differences,,  is  the  slight  directionality  associated 
with  the  acoustic  driver.  The  slight  offset  of  the 
laboratory  plot  nulls  is  caused  by  the  array  not  being 
exactly  aligned  to  broadside  at  0  Kohm  on  the  potentiometer. 
The  90  rotation  between  the  two  plots  is  caused  when  theta 
is  fixed  at  90  vice  0°  for  theoretical  plots.  Because  of 
the  narrowness  of  the  nulls,  laboratory  data  would  require  a 
sample  exactly  with  the  array  broadside  to  achieve  the  null. 
Other  than  the  difference  in  nulls,  there  is  exceptionally 
good  agreement  between  experimental  and  theoretical  beam 
patterns. 

In  the  investigation  of  three  element  arrays,  the  arrays 
were  configured  as  equilateral  triangles,  isosceles 
triangles  and  as  a  linear  array.  For  phase  shading,  three 
cases  were  looked  at:  all  elements  in  phase,  the  center 
element  180  out-  of-phase,  and  one  end  element  placed  180 
out-of- phase.  These  arrays  were  also  subjected  to  various 
amplitude  shading  (the  axial  quadruplet  composed  of  three 
elements  with  the  center  element  amplitude  shaded  to  twice 
the  value  of  each  of  the  end  elements  will  be  discussed  in 
the  quadruplet  section) ,  Experimental  data  for  a  three 
element  isosceles  triangle  with  all  elements  in  phase  and 
kDj^  =  kD2  =  2kD3  =  5.28  are  shown  in  Fig,  (4,7)  with  its 
theoretical  verification  shown  in  Fig,  (4,8),  The  deep 
theoretical  nulls  seen  on  Fig.  (4,8)  are  not  present  in  the 
laboratory  data.   This  can  be  attributed  to  a  small  variance 

45 


QQ 


UJ 


a. 


<: 

LU 
DO 

LU 


O 

CL 
I — I 

cr 


/ 


/ 


/ 


/ 


.■  / 


^.,:. 


/  . 


/ 


/" 


K-r-f- 


-te — R- 


\ 


•v .\:^ 

\  \  /N 

\  \ 

\ 


\ 


N 


S 


\ 


\ 


\ 


■v..      ^ 

)    I    IM 


.     /•• ../ 

•.         /•.  /  / 

/  / 

/ 


ro 


II 

CM 


/       S 


cn 
I 


m 

a 

a 


LU 

_i 
o 

•< 


o 
a 

> 

in 


CO 


a. 

•X. 
•< 


cn 
o 


a 
I 

UJ 

< 

X 

a. 


o 
cn 

0) 

rH 
0) 

u 

o 
cn 


-p 
c 


pa 

0) 


O 
-P 

u 


(U 


•H 
P4 


46 


THETA  VS.  DIRECTIVITY 

0 


PHI=  90  DEG3^iEE 

kDi-kD2=2kD3-  5.276 

D1=^3D2=0.3 
D3-0 15  I>4=0. 
P2=0  P3=0  P4=0 
A=l.  B-1.  C=0. 


Figure  4.8  Three  Elenent  Isosceles  Plot 
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in  the  equal  amplitude  shading,  a  small  amount  of  scatter 
from  the  large  array  center  joint,  and  a  slight 
directionality  of  the  individual  acoustic  drivers. 

The  beam  patterns  of  a  triplet  in  the  form  of  an 
isosceles  triangle  are  shown  in  Figs.  (4.9),  (4.10),  and 
(4.11)  to  show  the  effects  of  tilting  the  array  out  of  the 
horizontal  to  angles  of  theta  equal  to  30  ,  60  and  90 
respectively.   The  plane  of  the  elements  was   initially 

perpendicular  to  the  x-y  plane,   making   theta,   the  angle 

o 
between   the  z-axis  and  the  array  plane  equal  to  0    .   The 

array  plane  was  then  tilted  30  about  the  x-axis  and  rotated 

360   about  the  z-axis  to  obtain  the  data  for  Fig.   (4.9). 

The  array  plane  was  then  tilted  in  30   increments  to  obtain 

Figs.  (4.10)  and  (4.11).   For  this  triplet  case,  kD   =  kD ^  = 

2kD3  =  5.28  with  kD3  sensing  the  elevation  change  (Fig. 
(4.12)). 
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0.3 


0.3 


Figure  4.12  Array  Configuration  for  Tilt  Theta 


These  variations  in  tilt  angle  (theta)  show  that  the  3 
tilt  has  a  diminished  dB  level  for  all  lobes  except  the 
broad  lobe  coining  into  and  exiting  broadside.  This  is  due, 
in  large  part,  to  the  directionality  of  the  drivers,  because 
they  are  facing  the  microphone  to  a  greater  extent  at  3 
than  at  90  of  tilt  and  mutual  shading  exists. 
Additionally,  for  these  various  tilt  angles,  the  Z-axis  is 
rotating  about  the  local  vertical,  as  explained  in  section  D 
of  this  chapter. 

The  last  major  configuration  investigated  was  the 
quadruplet.  It  was  looked  at  both  as  an  axial  quadruplet. 
Fig.  (2.5),  and  as  a  planar  array,  Fig.  (2.7).  An 
experimental  beam  pattern  for  a  planar  array  quadruplet  with 


kD 


1   =  kD2  =  kD3  =  kD4  =  5.26  and  all  elements  in   phase   is 
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shown  in  Fig.  (4.13)  along  with  a  theoretical  beam  pattern, 
Fig.  (4.14).  There  is  good  agreement  between  the  two  plots. 
A  plot  of  an  axial  quadruplet  represented  by  an 
amplitude  shaded  triplet,  obtained  in  the  laboratory  is 
shown  in  Fig.  (4.15).  Comparing  Fig,  (4.15)  to  Fig.  (4.13), 
it  is  seen  that  when  two  elements  of  the  planar  array  are  in 
line,  the  arrival  from  the  more  distant  source  is  weakened 
by  the  shadowing  effect  of  the  nearer  source.  Therefore,  the 
near  source  can  be  emphasised  by  as  much  as  one  dB  when 
compared  to  the  far  source.  Additional  plots  of  various 
configurations  are  contained  in  Appendix  C,  where  good 
agreement  continues  to  be  seen. 

D.   EFFECTS  OF  ARRAY  TILT,  THETA 

Referring  back  to  Fig.  (4.2),  it  is  again  noted  that 
contour  lines  are  directly  related  to  polar  plots  when 
either  theta  or  psi  is  fixed.  However,  the  countour 
representation  is  not  as  useful  when  theta  and  psi  are  both 
changing  with  rotation  of  the  array  as  in  the  case  where  the 
array  has  been  tilted  in  theta.  In  this  case  the  polar 
curve  is  found  by  superimposing  a  trajectory  which  is  the 
specific  locus  of  0  and  ^  in  the  0  -"^  plane  onto  the 
contoured  surface. 

For  this  work,  the  microphone  was  fixed  in  the  far  field 
at  a  distance  of  3,2  meters  and  level  with  the  center  joint 
of  the  array.   Tilting  the  array  in  angle  theta  caused   an 
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THCTA  VS.  DIRECTIVITY 


PHI=  90  DEGREE 

kD^=kD2=kD3=l<D4=  5.275 

Dl=0 5  D2-05 
D3=0 5  D4=05 
P2=  0   P3=  0   P4=  0 
A=l.   B-0.  C=0. 


Figure  4.14  Theoretical  Planar  Quadruplet  Plot 
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offset  between  the  Z-axis  and  the  local  vertical  about  which 
the  array  rotated.  This  can  be  seen  in  Fig.  (4.16).  In  the 
left  sketch,  the  array  elements  are  aligned  in  the  standard 
rectangular  coordinate  system,  with  the  elements  in  the  x-y 
plane  and  rotating  about  the  z-axis.  However,  when  the 
plane  of  the  array  is  tilted  30^  about  the  x-axis,  Y'  =  Y+30 
and  z'  =  z+30°.  A  90°  rotation  about  the  z-axis  to  270 
yields  z  '  =  90°  -  30°  =  60°  .  This  means  that  a  36^ 
rotation  about  the  z-axis  causes  Yl  and  Zl  to  vary  from  120 
to  60  and  back  to  120  .  This  is  due  to  the  fact  that  the 
microphone  is  fixed  in  space  in  the  far-field  rather  than 
tilting  with  the  plane  of  the  array. 
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Figure  4.16   Z-axis  Offset  Effects 

This  problem  does  not  occur  in  the  classical  representation 
where  the  microphone  is  fixed  with  respect  to  the  axis  of 
rotation  of  the  array. 


E.   ERRORS 

The  sources  of  errors  in  the  beam  pattern  measuring 
system  were  examined  both  mathematically  and  empirically. 
Sources  of  error  for  direct  error  analysis  were  based  on 
multimeter  values  and  power  amplifier  settings.  With  output 
reading  from  the  multimeters  and  power  amplifier  exhibiting 
precision  to  the  third  decimal  place,  instrumentation  errors 
from  these  sources  are  less  than  1%  and  are  not  major 
contributing  errors  to  measured  beam  patterns. 
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The  major  source  of  error  in  the  system  was  determined 
empirically  to  be  the  acoustic  driver  positioning  on  the 
array.  The  effects  of  a  small  change  in  kd  are  shown  in 
Figs.  (4.17),  (4.18),  and  (4.19).  Here,  it  is  seen  that 
small  changes  in  kDj^  prove  to  have  a  large  effect  on  the 
size  of  the  small  lobe  at  ±  90°  and  0°  and  180°  ,  thus 
yielding  the  largest  error.  A  change  of  ±  5%  in  kd  effected 
the  structure  by  up  to  6  dB  as  seen  from  the  behavior  of  the 
small  lobe  at  i  90  .  This  change  relates  to  a  positioning 
error  of  about  ±  3  centimeters  and  presented  a  small  problem 
with  respect  to  exact  positioning  of  the  speakers. 

F.   VALIDATION  OF  RESULTS 

The  laboratory  beam  patterns  prove  to  be  in  good 
agreement  with  theoretical  beam  patterns.  Small  uniform 
beam  pattern  amplitude  differences  are  merely  a  consequence 
of  the  difference  in  normalization  methods  employed  in  the 
theoretical  and  experimental  beam  pattern  processing.  The 
theoretical  program  normalized  the  directivity  to  the  sum  of 
the  amplitude  weights,  whereas  the  laboratory  program 
normalized  the  beam  pattern  to  the  highest  received  voltage. 

An  independent  verification  of  the  experimental  program 
was  performed  by  a  thesis  group  under  the  auspices  of 
Professor  Don  Walters  of  the  Naval  Postgraduate  School 
Physics  Department,   who  asked  to  use  this   experimental 
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THETA  VS.  DIRECTIVITY 

0 


PH1=  60  DEGREE 

kD]^=kD2=kD3=kD4=  4.75 

Dl=0.5  D2=05 
D3=0.5  D4=0.5 
P2=  0  P3=  0  P4=  0 
A=l.  B=l.  C=l. 


Figure  4.17  kD^=  4.75  for  Error  Conputation 


60 


THETA  VS.  DIRECTIVITY 

0 


PHI=  60  EfiGREE 

U1=0  5  D2=0.5 
D3=0.5   D4=0.5 
P2=  0   P3=  0   P4=  0 
A=l.   B=l.  C=l. 


Figure  4.18  kD^^  5.0   for  Error  Conputation 
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THETA  VS.  DIRECTIVITY 

0 


PHI=  60  DEGREE 

kD^=kD2=kD3=kD4=  5.25 
D1=05D2=05 
D3=05   D4=0.5 
P2=  0  P3=  0   P4=  0 
A=l.   B=l.   C=l. 


Figure  4.19  kD2.=  5.25  for  Error  Conputation 
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system  for  related  work.  The  laboratory  system  was  placed 
in  use  by  this  research  group  in  the  selection  of  drivers 
for  a  5  X  5  array  currently  under  investigation  and  will  be 
further  utilized  to  establish  beam  patterns  for  this  array. 
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TABLE  II 
ARRAY  CONFIGnRATIONS 


Run  So/ 

Ampli 

tude 

Spacing 

No. 

Preq 

Phase    Tilt 

Shadl 

ng 

from  Center 

Elements 

:  (Hz) 

I  (deg)   1  (deg)   j 

(VAC) 

: 

(CM) 

1 

300 

a 

1 

300 

45 

1 

300 

90 

1 

9S0 

0 

1 

960 

45 

1 

960 

90 

2 

300 

0,0 

3 

30 

2 

300 

0,180 

0 

38 

2 

960 

0,0 

0 

30 

2 

960 

0,180 

3 

30 

2 

960 

0,0 

90 

30 

2 

960 

0,180 

90 

30 

2 

300 

0,0 

90 

30 

2 

300 

0,180 

90 

30 

2 

300 

0,0 

45 

30 

2 

300 

0,180 

45 

30 

2 

960 

0,0 

45 

30 

2 

960 

0,180 

45 

30 

4 

960 

0,0,0,0 

0 

30 

4 

300 

0,0,0,0 

0 

30 

4 

960 

0,180,0,180 

0 

30 

4 

300 

0,180,0,180 

0 

30 

4 

960 

0,0,0,0 

33 

30 

4 

300 

0,0,0,0 

30 

30 

Quadll/4 

300 

0,130,0, 180 

38 

30 

Quadl2/4 

300 

0,180,0,180 

60 

30 

Quadl3/4 

960 

0,180,0,180 

60 

30 

Quadl4/4 

960 

0,0,0,0 

60 

30 

Quadl5/4 

300 

0,0,0,0 

60 

30 

Quadl6/4 

960 

0,0,0,0 

90 

30 

Quadl7/4 

960 

0,180,0,130 

90 

33 

Quadl8/4 

960 

0,180,0,0 

0 

30,30,30,0 

Tril/3 

960 

0,0,0 

0 

30 

Tri2/3 

300 

0,0,0 

3 

30 

Tri3/3 

300 

0,180,0 

3 

30 

Trl4/3 

960 

0,180,0 

3 

30 

Tri5/3 

960 

0,0,180 

0 

30 

Tri6/3 

960 

0,3,180 

30 

30 

Tri7'3 

960 

0,180,0 

30 

33 

Tri8/3 

960 

0,0,0 

30 

30 

Tri9/3 

960 

0,0,0 

30 

30,15,30 

Triia/3 

960 

0,0,0 

0 

30.15,30 

Trlll/3 

960 

0,180,0 

0 

30,15,30 

Tril2/3 

960 

0,180,0 

30 

30,15,30 

Tril31/3 

960 

0,180,0 

60 

30,15,30 

Trll4/3 

960 

0,0,3 

60 

30,15,30 

Tril5/3 

300 

0,0,0 

60 

30,15,30 

Tril6/3 

960 

0,0,0 

90 

30,15,30 

Trll7/3 

960 

0,130,0 

90 

30,15,30 

Tril8/3 

960 

0,0,180 

90 

30,15,30 

Tri2a/3 

960 

0,0,0 

0 

,1, 

30,0,30 

Tri21/3 

960 

0,180,0 

0 

,1, 

30,0,30 

Tri22/3 

960 

0,0,3 

0 

,2, 

30,3,30 

Tri23/3 

960 

0,180,0 

3 

,2, 

30,0,30 

Tri24/3 

960 

0,0,130 

9 

,2, 

30,0,30 

Tri25/3 

960 

3,0,0 

0 

,2, 

30,30,30 

Tri26/3 

960 

0,0,0 

30 

-2, 

30,30,30 

Tri27/3 

960 

0,0,0 

60 

p2, 

33,30,30 

Tri28/3 

960 

0,0,0 

90 

,2, 

33,33,30 

Trl3a/3 

960 

0,0,0 

90 

>4, 

30,33,30 

Trl31/3 

960 

0,180,0 

90 

,4, 

33,33,33 

Trl32/3 

960 

0,3,3 

0 

,2, 

20,20,20 

Tri33/3 

960 

3,180,0 

0 

.2, 

20,20,20 

Tri34/3 

960 

3,130,3 

0 

-2. 

20,10,20 

Trl36/3 

960 

0,3,3 

3 

.2, 

20,13,20 

Tri37/3 

960 

0,3,3 

0 

.2, 

10,10,10 

Tri38/3 

960 

0,183,0 

0 

.2, 

10,10,10 

Tri39/3 

960 

0,130,0 

90 

.1. 

20,30,20 

Tri40/3 

960 

0,180,0 

90 

.2, 

.1      17, 

,9,17.9,17.9 

Tri41/3 

960 

3,0,0 

93 

.2, 

,1   17, 

,9,17.9,17.9 

Trl43/3 

960 

0,0,3 

90 

.2, 

,1   26, 

,8,25.8,26.8 

Trl44/3 

960 

3,180,0 

90 

,2. 

,1   26, 

.3,26.8,26.3 

Notet   (1)  Additional  data  runs  were  made  in  some 

of  the  same  configurations  with  fiberglass 
wrapped  around  the  centerpiece  to  identify 
the  effects  of  scattering. 

(2)  Additional  data  runs  were  mad«  In  some 
of  the  same  configurations  to  verify 
repeatadi lity. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  doublet  forms  the  basis  for  countless  sonar  and 
radar  arrays  currently  utilized  in  both  military  and 
civilian  applications.  It  is  the  ability  to  control  the 
amplitudes  and  phases  of  the  element  of  an  array  which 
allows  the  tailoring  of  the  array  radiation  pattern  to 
achieve  the  desired  element  of  interest  and  promotes  the 
maximization  of  array  performance  by  increasing  signal, 
decreasing  noise  or  both. 

The  ability  to  rapidly  measure  the  beam  pattern  of  an 
array  with  complete  confidence  in  the  results  is  the  basis 
behind  this  study.  A  laboratory  system  was  developed  to 
allow  collection  of  a  360  point  beam  pattern  in  a  45  to  50 
minute  period.  Software  is  easily  modified  to  allow  for 
collection  of  more  or  less  than  360  data  points.  One  of  the 
major  advantages  of  the  computer-controlled  method  of 
obtaining  a  beam  pattern  is  the  relatively  small  amount  of 
time  required  to  achieve  an  accurate  beam  pattern  based  on  a 
large  data  base. 

Recommendations  for  improving  the  current  laboratory 
system  are  based  on  reducing  electrical  noise  of  the 
rotation  motor  and  reducing  the  size  and  weight  of  the 
array.  First,  the  noise  of  the  rotation  motor  prohibits 
continuous   rotation  while   reading  the  resistance  on  the 
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potentiometer.  Since  a  CW  signal  does  not  require  broadband 
detection,  filtering  could  be  accomplished  by  use  of  a  lock- 
in,  fast  fourier  transform  or  heterodyne  filter.  Second,  the 
array  weight  should  be  reduced  to  lessen  the  drag  on  the 
rotation  motor  and  to  reduce  the  amount  of  coast-down  time 
required  for  heavier  arrays.  Next,  the  array  center  joint 
should  be  reduced  in  size  to  lessen  any  shadowing  effects 
from  the  large,  metal  center  joint.  (Scattering  effects 
seemed  to  be  small  judging  from  the  lack  of  any  effect  on 
the  polar  plots  when  this  joint  was  wrapped  in  fiberglass.) 
Also,  driver  mounts  need  to  be  designed  so  that  the  center 
of  the  driver  face  is  directly  on  the  axis.  The  2 
centimeter  offset  utilized  for  this  investigation  was 
minimized  by  making  all  readings  in  the  far  field.  Last, 
greater  distance  between  the  driver  and  microphone  would 
also  help  to  lessen  the  1  dB  variation  seen  when  an  array 
was  in  an  endfire  orientation.  This  only  affected  the  depth 
of  the  nulls. 
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APPENDIX  A 
EQUIPMENT  DESCRIPTION 

1.  ANECHOIC  CHAMBER 

All  experimentation  was  conducted  in  the  anechoic 
chamber,  Room  019,  Spanagel  Hall,  Naval  Postgraduate  School. 
The  anechoic  chamber  is  a  99%  echo  free  enclosure  for  sounds 
above  100  Hertz.  Acoustics  research  can  be  conducted  on 
sound  sources,  receivers  and  scatterers  with  minimal 
interference  from  surface  reflections  or  external  noises. 
The  size  of  the  outer  room  in  which  the  anechoic  chamber 
floats  is  37  feet  x  24  feet  x  21  feet  high.  Noise 
transmission  from  outside  is  minimized  by  the  floating  room 
within  a  room  construction.  The  outer  12  inch  concrete- 
walled  room  is  separated  from  the  inner  room  of  concrete- 
block  sides  and  floor  by  a,  2  inch  blanket  of  fiberglass  or 
cork.  Forty  inch  fiberglass  wedges  are  attached  to  all 
surfaces  of  the  inner  room.  These  are  designed  to  trap  and 
absorb  sounds  which  are  incident  upon  them.  A  suspended 
floor  is  constructed  of  a  grid  of  225  wire  cables,  each 
stretched  at  a  tension  of  150  to  200  pounds. 

2.  HP-86B  PERSONAL  COMPUTER 

The  HP-86B  computer,  operating  in  Hewlett-Packard  BASIC 
language,   is   equipped  with  a  64  Kilobyte   (K-byte)   memory 
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with  two  additional  128  K-byte  memory  modules.  The  system 
proved  more  than  sufficient  for  both  program  development  and 
as  an  instrument  controller.  Instrument  control  by  means  of 
the  computer  was  accomplished  with  the  HP-IB  (Hewlett- 
Packard  Interface  Bus)  which  is  a  standard  connection  on 
most  HP  equipment.  This  facilitated  instrument  control  for 
the  various  system  components. 

3.  HP  7470A  GRAPHICS  PLOTTER 

The  HP  7470A  Graphics  Plotter  is  a  vector  plotter  which 
utilizes  HP's  micro-grip  drive  technology,  moving  both  the 
paper  and  the  pen.  Multicolor  plotting  is  available,  but 
there  are  only  two  pen  positions  so  pauses  must  be  used  to 
change  pen  colors.  Only  one  line  of  program  code  is 
required  to  shift  the  output  to  and  from  the  plotter  to  the 
HP-86B  CRT. 

4.  HP-4192A  LOW  FREQUENCY  IMPEDANCE  ANALYZER 

The  HP-4192A  LF  Impedance  Analyzer  is  a  fully  automated 
instrument  designed  to  measure  a  large  range  of  impedance 
related  characteristics.  It  can  be  used  to  measure  eleven 
impedance  parameters:  absolute  value  of  impedance  Z, 
absolute  value  of  admittance  Y,  phase  angle Gf  resistance  R, 
reactance  X,  conductance  G,  susceptance  B,  inductance  L, 
capacitance  C,  dissipation  factor  D,  and  quality  factor  Q. 
The   HP-4192A  was  used  to  measure  the   frequence   response. 
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impedance,  resistance  and  reactance  of  the  University  Sound 
ID-30  driver  units  as  discussed  in  section  10  of  this 
appendix . 

5.  HP-3421A  DATA  ACQUISITION/CONTROL  UNIT 

The  HP-3421A  Data  Acquisition/Control  Unit  (DAC)  has  a 
5.5  digit  voltmeter,  a  10  kHz  counter  and  an  HP-IB 
controller  interface.  The  DAC  was  utilized  as  a  computer 
controlled  switching  unit  with  channels  1  and  2  configured 
as  actuators.  This  provided  a  voltage  path  to  control 
relays  in  the  motor  power  and  direction  controller  for 
rotation  of  the  ID-30  speaker  array,  Fig.  (A.l) . 

6.  HP-3478A  MULTIMETER 

The  HP-3478A  Multimeter  provides  both  resistance  and 
voltage  measurement  capabilities  with  3.5  to  5.5  digit 
resolution  with  +  0.3  mV  accuracy  on  the  100  Hz  to  20  kHz 
scale.  The  multimeter  was  remotely  controlled  by  the 
computer  via  the  HP-IB.  Voltage  readings  from  the  GR-1962 
microphone  and  resistance  readings  from  the  variable 
potentiometer  on  the  array  rotator  being  read  by  the 
multimeter  and  then  sent  to  the  HP-86B  computer  via  the  HP- 
IB. 

7.  HP-3314A  FUNCTION  GENERATOR 

The  HP-3314A  Function  Generator  is  a  multimode  function 
generator   which  provides   sine,    square,   and   triangular 
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Data  Acquisition  /  Control  Unit 


I 


Q 


1         23456789 
QOOOOOOOO 

+  5  Volts 


© 

+  5  Volt 

Power  Supply 

0 


oooooooo 
Common 


6 
Ch 


to  HP  Multimeter 
Ohm  sensing 


6 
Ch  2 

O 


^ 


to  Motor  Powers 
Direction  Controller 


Figure  A.l.   Data  acquisition  connection  to  motor-power 
controller. 
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functions  as  well  as  any  desired  arbitrary  wave  form  from 
frequencies  of  0.001  Hz  to  19.999  MHz.  It  was  operated  in 
the  free  run  mode  providing  a  continuous  1  Vac  sine  wave  to 
Hp-467A  Amplifier's  inputs. 

8.  HP-465A  AMPLIFIER 

The  HP-465A  Amplifier  is  a  general  purpose  preamplifier 
and  impedance  converter  from  10  megohms  to  50  ohms.  It  has 
a  selectable  gain  of  20  or  40  decibels  and  is  stable  over  a 
frequency  range  of  5  Hz  to  1  MHz.  The  amplifier  was  used  to 
amplify  the  received  signal  from  the  GR-1962  microphone. 

9.  HP-467A  POWER  AMPLIFIER 

The  HP-467A  Power  Amplifier  is  a  10  watt  peak  power 

amplifier  and   a  -20  to  +20  void  do  power  supply.    It  has  a 

wide  band  width  with  low  dc  drift  from  5  Hz  to  1  MHz.  A 
variable  gain  control  allows  adjustments  of  gain  between 
zero  and  10  with  resolution  better  than  0,1%  of  full  power. 
Four  HP-467A's,  one  for  each  ID-30  driver  unit,  were 
utilized  in  adjusting  power  inputs  from  the  function 
generator . 

10.  UNIVERSITY  SOUND  ID-30  DRIVER  UNITS 

The  University  Sound  ID-30's  are  medium  power  horn 
driver  units  designed  specifically  for  wide  midrange 
response.  Their  tropicalized  voice  coils  are  2  inches  in 
diameter  with  linen  base  molded  phenolic  diaphragms  with  a 
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minimal   16  ohm  impedance.    The  frequency  response  is  85  Hz 
to  7.5  Khz.   They  were  used  as  the  sound  sources. 

11.  GENERAL  RADIO  TYPE  1962  MICROPHONE 

The  GR-1962  microphone  is  a  0.5  inch  diameter  electret 
condenser,  microphone,  providing  a  flat  response  at  zero 
decibels  between  20  Hz  and  15  Khz.  It  experiences  rapid 
drop  off  above  15  Khz  and  is  nonresponsive  above  21  Khz.  A 
General  Radio  Model  1560-P62  inline  power  supply  was 
connected  to  the  microphone  1560-P4J-  preamplifier  to 
amplify  the  received  signal. 

12.  HP-6215A  POWER  SUPPLY 

The  HP-6215A  power  supply  provides  zero  to  30  volts  dc 
or  0  to  500  Mamps .  It  was  used  to  provide  switching  voltage 
for  the  Motor  power  relay  in  the  direction  controller. 

13.  MOTOR  POWER  AND  DIRECTION  CONTROLLER 

The  controller  for  the  array  rotator  controlled  power  to 
the  rotator  and  also  the  direction  of  rotation.  With  leads 
to  the  rotator  assembly  potentiometer,  it  also  provided 
resistances  which  corresponded  to  a  given  angle  (0  to  10 
Kohms  translated  into  0  to  360  degrees) .  A  +5  volt  dc 
signal  to  IC-1  via  channel  1  of  the  DAC  provided  on/off 
power  to  the  rotator  and  +5  Volts  to  IC-2  via  channel  zero 
of  the  DAC  provided  clockwise  rotation  of  the  driver 
assembly  unit  (Figure  A. 2) .  [Ref .  8] 
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Dip  R 

elay  Input  and  Output 

1. 

AC  high 

2. 

AC  low 

3. 

+5  V  from  DAC  Chi 

4. 

Common 

5. 

+5  V  from  DAC  Ch  0 

6. 

Ohm  sensing  High 

7. 

Ohm  sensing  Low 

Power  to  Motor 


a)  Circuit  Diagram  for  Modified  Motor  Power  and  Direction  Controller 


b)   360"  /  10  kO  Potentiometer  c)    16  Pin  DIP  Relay 


Figure  A. 2.   Circuit  diagram  for  modified  motor-power/  direction 
controller  [Ref.  9.,  p.  57]. 
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APPEM)IX  B 

BEAM  PATTERN  COMPUTATION  AND  PLOT  PPDGRAM 

10  pr![MT  " yhnrl  (■mFTLRN  COMf-UIATION  r.MU  FM-uT  r'RuGPAM  .♦♦»♦" 

1  1  FRIfJI  "«•«»•»»•«♦♦»«••♦■»♦♦♦»♦•#•♦«♦••»♦♦♦«»♦»♦««•#*♦»«•♦••«" 

1 :   ' 

13  I 

Z'3  PRINT  ".*••♦  SET_UP  »»•♦•" 

22  • 

30  PR  IMF  CHR$  \  H3  ) 

40  PRIliT  "rillS  PROGRAM  CONTROLS  THE  FOLLOWING  APPARATUS  TO  MEASURE,  DISPLAY," 

53  PRlfil  "PLOl  hND  STORE  THE  DIRECTIONAL  CHARACTERISTICS  OF  A  SINGLET," 

&0  PR  INF  "DOUBLET,  TRIPLET,  AND  OUnDPUPLET  ARRAY." 

Bl  I 

e:    I 

70  PRINT  "  INSTRUMENT  ADDRESS  " 

00  PR  INI  •' " 

93  PRINT  "  HP.  3473  MULTIMETER  ( UOLTS  )     724 

100  f'RINI  "  HP_347B  MULTIMEIER  '.RES)       723 

no  PRINT  "  DnJh    AQU  1  S  I T  I  ON  '  CONTROL       709 

111  I 

12^  DiSP  "INITIALIZE  INSTRUMENFS" 

121  I 

130  REM  •»♦  PRESET  AND  INITIALIZE  INSTRUMENTS* * • 

140  OiiTPUl  709  ;  "RS"  i  RESET 

150  OUTPUF  723  ;"r3R5Tl"  i  FREE  RUN  30K0HM5 

lEO  uUTf'UT  724  ;"F2RATt"  •  FREE  RUN  AUTOPr.NGE  AC  UOLTS 

151  I 

170  DISP  "INITIALIZE  PROGRAM" 

1711 

180    OPTIOFJ    EASE    1 

i '00    CiJM    nliCLEUbO  )  ,CALC_ANGLE(  150  )  .NAME*[30] 

200    DEu 

2  10  PRINTER  IS  701 
220  INFEOfR  I ,C ,L ,N 

230  Dirl  PJ  360  '<  ,UOLT(  ZLQ  )  ,OHMS_DATA(  3S0  )  ,X(  3G0  )  ,Yv  360  )  ,K\  360  ) 

2  40  M"uL[-0 

250  MASS  STORAGE  IS  ".DATA  I" 

260  ni?P  "NUMBER  OF  DATA  POINTS  DESIRED?" 

2/0    BEEP    10,300 

2&0    I Nf 'UT    M 

2  90  DtSP  "UANI  NEU  DATA  (Y/N)"?" 

300  BEEP  10  'SOd 
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Z\v}    INPUT  0$ 

12^'    11  0*-"  r  "  THEN  bOTO  550 

.■:.:0  IF  i>l  =  "M"  TMFN  GOTO  700 

340  GOTO  300 

3^30    f'Eli    ►•»    OEFnULT    PnRr>,(-lhTEF5    •♦» 

33  1     I 

333     I    hOsT    INSERT    FuR    LINES    390,    400,    410    mNU    610    PRIOR     fO    SThRT    OF    RUN 


3'30 
4  00 
410 
4  30 
4::0 
440 
4  SO 
4'.-)0 
^&^^ 
fi  1  0 

53  0 
530 

54  0 
550 
5b  0 
570 
575 
5b0 
5 '30 
600 
BIO 
GZO 
G30 
G40 
550 
Ei^O 

b-:)0 

700 
■/ 1  0 

73  0 
730 

74  0 
750 
7G0 


CREfiTE    ■■♦    INSERT    NftME    OF    RUN    •",350.16 

AbSIGNtt     1     TO       "»    INSERT    NriNE    OF    RUN    •" 

Nr,riE$-"»     iriSERT    NhME    OF    FUE    *"     i     FILE    NAME 

FULL  SrMLE  =  500    i    L0CK_IN    f-VJnC 

l0(IU  =  3G0/ 10000    I    lONUERT    OUflS    TO    DEGREES 

SET__UP:     I     *•«    5ET_.UP    MErVGUREMENT    PnRr.HETERS    »•• 

FRiNf    GHR$    U0  ) 

DIGF  "MhNUALLY  ENTER  DESIREO  FREQUENCY  FOR  SPEAKER  DRIUER" 

DI3P  ■■ENTER  O  IF  NO  CHANGE" 

DI3P  "ENFER  I  TO  GHANT^E  1  ITL  E" 

EEEP  10,600 

INPUr  G  !  CHOICE 

IF  G=0  THEN  GOTO  560 

IF  G-1  I  HEN  GOTO  750 

PRINT  GHRi  (  lO  ) 

OUTPUT     ,'09    ;"GL50"     i     ROTATE    CW 

RRINI     "  I"  ,  "OHflb^DATA    ujHN^  )  '  ,  "  UOlT    (  UAC  r' 

FOR  1=1    ro  N 

OUTPUT    703    i"CL3r'     i     POWER    TO    MOTOR    ON 

WAI]     .5     I    MILLISECONDS 

OUTPUT    709    ;"0PN1"     i     P'llJER    TO    MOTOR    OFF 

UMlT    4  000 

ENTER    73  3    ;    OHMS_.DATA<,  I  ) 

LfllER    ','34    ;    UOLT(  I  ) 

PRIfir     I  ,OHMS  _DmTA(  I  )  ,UOLT(  I  ) 

NExl     I 

PGP     f  =  l     Tt)    N 

PRlNTtt     1     ;    OHMS_DATA(  I  )  ,UOL  T(  I  ) 

NE'',  r  [ 

ASSIGNfl     1    TO       • 

GOSUB  PLUT_.ROUTINE 

GOTO  7'/0 

ENlEf<    NfiMES     I    MEN    TITLEi 

bul  u    4ijO 
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770  LNlJ 

■I'J^  ri[5F    "WHnr    15    DATA    RUN    DESIRED'^" 

790  BEEf^    10  ,500 

S\-d  rtf.5IG(ltt     !     TO       "•     iriSERT    rinhE    OF    RUM    ♦" 

S.:0  FOR    I- 1    TO    fj 

03iZ)  RLriOit     I     ;    UHMS_DATPi(  I  )  ,UOLT(  I  ) 

5-iO  rJE/l     I 

B50  Ml:5I6lltt    1     rO       » 

3b 0  GOTO    -^z.o 

861  I 

aG2  I 

B70  PLOr_ROUTINE:     i     ♦»•    PLOTS    FOLnR    DPiTn    »♦* 

B7  1  I 

8  7:   I 

8B0  PLOTTER  15  705 

C'33    GRnPH 

900    Gll£,''iR 

91 '2)    PEli    1 

920    BCmLE      R0  ,60,-60,60 

930    XfiAlS    0,10  ,  -40  ,40 

94  0    YfiMG    0,10,-40,40 

950    CGIZE    3 

960   LltJL    ]  rPE    1 

970  FOR  L=0  TO  40  STEP  10 

960  r-hi'.^E  0  ,L 

990  1  >')iv'j    b 

1000  L.-iPEl.  L-40 

1010  HE  XT  L 

1020  1  iriL  TYPE  1 

1030  C5IZE  6 

104  0  f-M"E  0  ,50 

10  4  1  ! 

1012  I  HLiERT  LABELS  1050,  1080  AND  1100  PRIOR  TO  START  OF  RUN 

1043  I 

1050  LABEL  ;"OUADRUPOLE  BEAM  PATTERN  (dB)" 

)0C0  GSIZE  3 

1070  li'V'E  0  ,  •  45 

UlciO    LnBEL    ;"R-3.20    M","  "  , "  I  d-3  .  30"  ,  "  ","FREQ  =  30O    Hr" 

1090    Mf.i"E    0,-50 

1100  LABEL  i  "PHASE  =  0,  180  ,0  ,  180  Deq","   ","AMPL=  0.093  Vac\" 

1110  Ni'i'L  45,0 

1120  Li;t-G  5 


."AI1GLE  =  30  Deg' 
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130    LABEL     ; "MIKE" 

1  4 'J    LiriE    TVF^E    4 

IE0    FOR    L=  115    TO    40    STEP     10 

IbO    |-|iii^E    I   ,0 

170  FOR  1=0  TO  3G0  STEP  5 

l*j0  DRAW  L'L.uS  n  )  ,L»Sirj  (  I  ) 

190  riEXT  I 

200  NEXT  L 

2  10  LIME  TYPE  I 
220  FEU  UP 

2  30  FEN  2 

240  FOR  1=1  TO  U    i  BEGINS  BEhN  PnTTERN  COMPUTnTIuN 

250  IF  OOLT^  1  ,■  MOuLT  THEfi  N'JOL  T  =  UOLT  ^  I  ) 

2G0  NEXT  I 

270  FOR  I-l  TO  N 

280  K=-uHliS_nhThi  I  )/lO000'3b0 

2  90  K 1=K 

300    IF    OflMS..  OATA(  I  )     10    ItlEN    K=K  1 -(■360/295 

310    IF    uHli5_DhIfin  )    10500    THEN    GOTO    1400 

12'^    '.'EC  lOh  =  20«LGT    C.'OL  J i  I  )/r-lWOLT  )  ^40 

330    LINE    TVFF    3 

340    IF    'JFlTOR    1     niEN    i,'ECTOR  =  0 

350    PRINT     I  ,K  ,'JEOTOR 

3G0    >0:  I  )-iJECTOR'COS    (  K  ) 

370    V(  I  i--l'ErTOR«SIN    (  K  ) 

SB\^)    PLOT    '<  (  1  )  ,Y(  I  ) 

390    PErj    UP 

400    NEXT    I 

4  10    RETURN 
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TABLE  III 
SAMPLE  PROGRAM  DATA  COLLECTION  OUTPUT 


4 

5 

6 

7 

8 

9 

10 

1  1 

12 

13 

14 

15 

IG 

17 

18 

19 

20 

21 
-)  -1 

T  "7 

24 
25 
2B 
27 
2S 
29 
30 
31 
32 
33 
34 
35 


OHPIS^DmTA    (0HM=  )              [>0i 

.T    (  UrtC  1 

4 

;  4 : :  4 

23 

'8029 

3G 

'5352 

56 

'02S3 

79                                                .F 

.6387 

106                                              .£ 

J  026 

12  3                                              .i 

w753 

160                                              .^ 

d4058 

194                                          '    .E 

d2055 

215                                              .E 

;i957 

251                                              .E 

d3613 

281                                              .i 

J5734 

ZiJti                                              .E 

]0736 

330                                              .E 

J5774 

359 

'2735 

401                                              .i 

i0624 

42  9                                              .i 

^6951 

458                                       .^: 

33089 

450                                                1 

00539 

522                                              1 

07  14  9 

543                                              1 

1  2759 

572                                              1 

155-77 

600                                                       1 

24  4 

628                                              1 

28525 

6-18                                              1 

'^.  "^  "5 

E73                                              1 

35  172 

704                                              1 

40042 

725                                              1 

42822 

755                                             1 

4  55  4 

786                                              1 

47833 

514                                              1 

43626 

646                                              1 

50753 

069                                             1 

51  197 

895                                              1 

5  11  00 

92  4                                              1 

50432 

78 


I 

T 


10 


15 
lb 
17 
13 

1  J 
20 

2  1 


2b 

2  3 
2  3 
30 
31 


TABLE  IV 
SAMPLE  CALCULATED  BEAM  PATTERN  OUTPUT 


;4 

^5  33.2G4 


Degrees 

1  .3b 4 33538305 

Decibels 

3  1  .  J  1  7.j537e27 

.d2  3 

34. 

25425-O403 

1  .236 

T  "7 

,5544395323 

2.01b 

J;  J  . 

3461423354 

2.844 

-!/  u.  . 

,3507973377 

3.81b 

-7  n 

232559503b 

4.E03 

31  , 

.6423331572 

5.7G 

31  . 

0727507067 

5.934 

30. 

,733552392 

^.74 

30. 

7220173553 

3.03G 

30. 

,9945387053 

10.116 

3  1  . 

4360046043 

1  1  .033 

32. 

.0852078323 

1  1  .33 

Jj  *.  . 

77022  1557  1 

13.234 

-1'  -J  . 

,6440053809 

14. 4  36 

34. 

5334233377 

15.444 

35. 

.  I94G23  1359 

16.438 

35. 

737 1 033  !  67 

17.L-G5 

3b  , 

.4553277425 

15.792 

-7-7 

,  00339896 1 3 

1  9  .  b  4  i; 

37  , 

.4521602313 

20.592 

~7 

.9133967977 

2  1.6 

3  3  , 

.305544 1246 

22.603 

33  . 

.609  1395402 

23.323 

38, 

.  837050 1316 

24.223 

39. 

.09039284 

25.344 

39 

.3343026072 

25. 136 

39, 

.5050337255 

27.  1  3 

39 

.6^47499227 

2B.29b 

39, 

.3045643216 

2  9.30  4 

39 

.9092773335 

30.456 

39, 

.9^47433642 

31  .234 

40 

T  ~   ~i  n 

39 , 

.  994355667 1 
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APPENDIX  C 


Experimental  and  Theoretical  Plots 
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THETA  VS.  DlRECllVITY 

0 


PH1=  90  DEGREE 

kDi=kD2=2kD3=-  5.276 

D1=K).3  D^=0.3 
D3-0 15  D-4=0. 
P2=100P3=0  P4=0 
A=l.  B=l.  OO. 


)(rEco 


Figure  C.2  Theoretical  Three  Element  Isosceles  Plot 
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THETA  VS.  DIRECTIVITY 

0 


0(tECO 


PHI=0  DEGREE 

kDi=kD2=kD3-kD4^  5.276 

D1=0  5  D3=0  5 
D3=05  D4=05 

P2-0  P3=180P4=180 
A=l.  B=l.  C=i. 


Figure  C.4  Theoretical  Planar  Quadrupole  Plot 
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